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Abstract 

We propose here an alternative interpretation of the fermi-Unarization approach 
to interacting electron systems, based on the requirement that the coefRcients of the 
linearized operators are Clifford-like variables, whose anticommutator equals an un- 
known constant c. We apply the approximation to the Falicov-Kimball model, explic- 
itly solving the self-consistency equation for the unknown, which turns out to behave 
as an order parameter. We discuss its relation with a metal-insulator transition and 
some thermodynamical quantities. In particular we show that our approximation in 
the T = limit reproduces exactly the Gutzwiller results for the Hubbard model. 
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1. Introduction 

Both systems of itinerant interacting electrons on an infinite lattice, and 
simple systems of few electrons interacting with bosons, are generally described 
by hamiltonians whose dynamical algebra is infinite dimensional. Various ap- 
proximation techniques have been developed in order to deal with such systems. 
In particular, in a set of recent papers^^"^' an approximation scheme was pro- 
posed, referred to as fermionic linearization scheme, which can be applied to a 
generic many-fermion hamiltonian. It consists in replacing in the hamiltonian 
certain bilinear products of (sums of) electron creation or annihilation opera- 
tors, say ai and 02, by terms linear in some fermion operator / multiplied by 
appropriate Grassmann-like coefficients 6, 

aia2 + ala\r^ef + fe , (1) 

where {6,9} = 0, and {9 J} = /} = - {9,ai} = {^,aj} {i = 1,2). The 
anticommutation relations of / and are uniquely determined by ai and a2, 
and depend on the problem studied. The fact that both the operators /, and 
the Grassmann coefficients 9, 6 satisfy anticommutation relations guarantees 
that the bilinear products on the r.h.s. of (1) have the same 'statistics' of the 
bilinear operators at the l.h.s.. 

Once substitution (1) is performed, the scheme allows one to obtain the 
spectrum of the linearized hamiltonian - after recognizing that the 'effective' 
model has a dynamical algebra which is a 2i2-graded algebra - via an inner 
automorphism of the algebra itself (which generalizes the customary Bogolubov 
rotation). 

In the present note, we propose a new view of the fermionic linearization 
scheme, which consists in requiring that the variables 9, 6 satisfy a Clifford-like 
instead of the Grassman-like algebra. More precisely, we set 

{9,e} = c^ , c e H , (2) 
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with c an undeterminate, to be defined for each specific problem. Notice that the 
requirement (2) on the ^'s imphes that the dynamical algebra of the linearized 
model is no longer graded, but simply a Lie algebra. In other words, we require 
that the ^'s behave as operators rather than anticommuting numbers. Indeed, 
by inspection of (1) one can easily verify that, in the simple case in which ai 
and 02 are single electron operators, and {/, /^} = 1, cqn. (1), with c = 1, 
maps a two-electrons operator into another two-electrons operator, hence the 
approximation of the r.h.s. term of (1) becomes exact. In general, this is 
not true, and a value of c has to be determined self-consistently according to 
eq. (1). The self-consistency equations reconduct then the exact results for 
the linearized model to approximate (mean-field like) results for the original 
hamiltonian. 

In the following, we will use this approximation to the solution of the 
Falicov- Kimball model. The latter gives a very simplified description of a system 
of itinerant fermions interacting only locally. In this case, prescription (1) is 
applied to the itinerant part of the hamiltonian, reducing it to an effective single- 
site operator, while it leaves unchanged the interaction term. The resulting 
approximation thus in principle goes well beyond the standard weak-coupling 
mean-field theory, and indeed it turns out to be capable of describing a metal- 
insulator transition. 

Let us observe that the approximation (2) was already used in a different 
context t"^] with a fixed value for c, i.e. c — 1. 

2. The Falicov-Kimball model 

The Falicov-Kimball model'^^ provides a very simple description of large 
systems of itinerant interacting fermions, by considering two different species 
of electrons (say with up and down spin) on a lattice A, one of each itinerates 
on A, the electrons with opposite spins being fixed at their sites, and assuming 
that the electrons interact only via an on-site Coulomb repulsion term. The 
grand-canonical hamiltonian reads 

^FK = -Mn$^iVi-/^dJ]A-2t J] ^Aj + C/J^iViA , (3) 

i i <i,j> i 
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where t > is the hopping amphtude, and J7 > is the local electron-electron 
repulsion. are operators which create and annihilate the itinerant elec- 

trons {{Ai,Ai} = , {AlAj} = Ni = AlAi, i, j e A), and < i, j > 

stands for non-oriented nearest neighbours (n.n.) in A. Moreover Di is the 
number operator of the non-itinerant electrons. As the operators Ni and 
Di both commute with the hamiltonian, the chemical potentials /in and fXd 
allow to fix the average number of electrons of the two species. 

The Falicov-Kimball model was introduced for studying the metal-insulator 
transition in transition metal oxides, and can be considered as a simplified 
version of the Hubbard modelt^L The exact statistical mechanical solution for 
the model described by Hpx is known only for large dimensions'^]. However, a 
few general theorems are knownl^l for the symmetric (or neutral) case /^n = A*d = 
— , and in particular an Ising-like phase transition is expected for dimension 
D > 2 at some critical temperature, whose value should vanish both for small 
and large U. Moreover, there are a number of investigations of the ground 
state phase diagram in dependence on the configuration of fixed spins'^!. Also, 
a strong-coupling {U » t) thermodynamic mean- field theory - based on the 
D = oo exact solution - was proposed'^^l . 

The fermionic linearization approach'^"^!, mentioned in the introduction, 
provides as well a powerful approximation scheme for Hubbard-like models in 
the strong coupling limit. In fact it treats in an exact way the Coulomb inter- 
action term, whereas it acts only on the hopping term. Let us write the latter 
as 

j:AlA, = lJ:{elA + Ale-^ , (4) 

<i,j> i 

with ©i = - \^ Ai, q denoting the number of nearest neighbours of a site 
q ^-^ 

jn.n.i 

in A. Of course, the operators Oi have non trivial anticommutation relations 
among themselves as well as with the Ai's. On the other hand, in [1-3] the 
©i's were approximated by variables ^i's anticommuting among themselves as 
well as with the fermion operators, i.e. {^i,^j} = = = 
= {^i, Aj},Vi,j e A. The former prescription is exact only for i and j far 
enough, depending on the lattice A, whereas it is definitely too simple for i and 
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j coinciding or having nearest neighbours in common. 

Here we propose to improve the fermionic approximation scheme by re- 
placing the operators ©i by variables still anticommuting with the fermion 
operators, and satisfying the following algebra (which is a straightforward gen- 
eralization of (2)) 

{0i,0j}^c^S,,i , {^i,^j} = . (5) 

Once the above approximation is inserted in (3), one obtains a reduced hamil- 
tonian Hfk which is a sum over lattice sites of single-particle hamiltonians, 7ii, 
commuting with each other, 

Hi = -iinNi - iidDi - tqc{f]iAi + Alrii) + UN^D;, , (6) 

with ?7i = — ,so that {f]i, = 5\,y 

The Di's arc to be considered as classical, Ising-like, variables, whose two 
possible eigenvalues and 1 label two orthogonal projections of Hi == 'hif'^ ® 
l-6^\ The problem of finding the spectrum of hamiltonian (3) is thus reduced, 
after linearization, to that of diagonalizing the local effective hamiltonian Ti.^^'-*. 
In order to do it, one should first identify the dynamical algebra, of (6); it 
is easily verified that the latter coincides with tt(2), generated by 

^ = «(2) = {iVi±f7ir7i;r7iA±4r/i} . (7) 

The transformation which rotates the hamiltonian into its diagonal form TIfk is 
then obtained by acting on 'h[^'^ with exp (adZ) = ^ ~['^' [^■' ■ ■ ■■> •] ■•■]]) 

n=0 

where Z is an appropriate skew-hermitian non-Cartan element of A, Z = 

t 2r 
p{r)iAi — Aj?7i). It is easily verified that the choice p = arctg— — , with 

(J Ui fiji 

T = cqt, implies 

Hfk = \ ^-Amm + A^i) ± \j ^\ ^ ^T'^{n,7]^ - - /^dA , (8) 
with e\ = \JD\ — /In', 'Hfk is manifestly diagonal. 
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The result (8) is also interesting from the point of view of statistical me- 
chanics, in that the partition function Z is immediately obtained from (8) as 

Z= exp(^-pnFK) . (9) 

Ni,Di,fjirii=0,l 

Predictions for physical quantities can then be obtained from Z once the average 
numbers of electrons of the two species are fixed through the chemical potentials, 
according to 

1 dZ 

n=<Ni>=—-— ■ (10.1) 

I3Z dfin 

I dZ 

d=<D;>=——- , (10.2) 

PZ dfid 

where < • > stays for the thermodynamical average in the Gibbs ensem- 
ble of operator • (i.e. < • >= Z!iVi,Di,^j?=o,i • ^^P = 
EjVi,A,w=o,i exp (adZ)(«) exp {-PHfk) )• 

Moreover, in order to have quantitative predictions, a numerical value for 
c has still to be self-consistently determined. Indeed, the prescription of substi- 
tuting in the hopping term the ©i operators with the ^i's can be implemented 
once more in (4) , giving rise to the following self-consistency equation, 

< fjAi + AIt) >= 2c < fjr] > , (11) 

in which we have assumed translational invariance of the lattice, implying rji = 
?7, Vi e A. 

The three equations (10.1), (10.2), and (11) have interesting features. First 
of all, we notice that (10.2) can be solved explicitly for /x^, and gives the result 

^ 1 + e^^" + 2e^ cosh ^ ^//^ + 4t2 

expf3iid = z ^ • (12) 

1-^1 + e^(M.-t/) + 2e^^^ cosh§V(/^,-C/)2 + 4r2 

Moreover, it is easy to check that eq. (11) always factorizes a solution c = 0, 
which correspond to the insulating behavior. Besides this solution, in general 
the system formed by (lO.l)-(ll), with //^ given by (12), is highly non-linear, 
and must be dealt with numerically. It turns out that it has different non-zero 
solutions. The physical one is to be chosen as that which minimizes the Gibbs 
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frc^energy /, / = -i In Z. In the next section we shall discuss the results of 

the numerical analysis, as well as the analytical results which can be obtained 
in some limiting cases. 

3. Results and discussion 

kT 

In figure 1. we report the mean- field parameter c vs. temperature — , at 

qt 

half- filling and for the symmetric case n = d= ^. In this case it is easy to check 
that the solution to (10.1)-(10.2) is iJ,n = l^d=^r- c is plotted for different U 
values, and exhibits a typical order-parameter like behavior. For U = (non- 
interacting case) it rises from zero, in the high-temperature regime, to one, at 
T = 0. For generic U < Aqt, it is possible to show rigorously that, in the limit 
T — > 0, c reaches a value cq given by 

where U = This suggests that the value c = 1 used in [4] is correct at half 
filling, only in a low-temperature non-interacting regime or for D = oo. On the 
contrary, for U > Aqt, the only solution to (11) is c = 0. 

The expression (13) for cq clarifies the physical meaning of the parameter 
c. Indeed, recalling that on a hypercubic lattice q is twice the dimension of the 
lattice, eq. (13) reproduces exactly the Gutzwiller result I-'^-'^] for the discontinuity 
in the single particle occupation number at the Fermi surface, obtained for the 
conventional Hubbard model when T = 0. This is not surprising as, on the 
one hand, the Gutzwiller result for the Hubbard model was obtained in fact by 
neglecting the kinetic energy of one species of electron, thus in an approximation 
very similar to that at the basis of the Falicov-Kimball model. On the other 
hand, according to eqs. (4), (5), and (11), at half filling c coincides with the 
expectation value of the hopping term, and hence is related to the discontinuity 
in its Fourier transform. 

Notice that when U = then cq = 1, and the ground state has all the 
electrons below the Fermi level. For any c 7^ 0, the ground state has some 
electrons above the Fermi level, but the gap is still there, and, according to eq. 
(5), the generic lattice site on which one has confined the linearized hamiltonian 
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is still interchanging fermions with the rest of the lattice. When cq = on the 
other hand, the gap in the density of states disappears, and at half-filling we 
have exactly one electron per site. In this case, the remaining of the lattice 
behaves as a system of scorrelated 'average' fermions (i.e. as if they were frozen 
at their own sites) and we are in presence of an insulating phase. 

The above analysis suggests that c could be able to describe the transition 
from a conducting to an insulating state. Indeed, again in agreement with the 
Gutzwiller result, at T = we find that the double occupancy expectation value, 
V =< NiDi >, vanishes precisely aX U = 4. Explicitly, analytic calculation 
shows that 

7^=14^^ ---- . (14) 

^ otherwise 
It is worth noticing that the result (14) coincides with the exact result both in 
the limit U = and in the limit U » 1. 

A deeper analysis of figure 1. shows that the transition from non-zero 
to vanishing c is of different order depending on the value of U. Indeed, by 
requiring that (11) vanishes also around |c| = 0, one can verify that there exists 
a tri-critical point at [/ — Ut, where Ut is solution of 

tanh — = y . (15) 

One finds a numerical value Ut — 1.845. For U smaller than Ut the transition is 
second order, and the critical temperature is found analytically as the solution 
Tc of the following equation (obtained by requiring that (11), upon factorizing 
the c = solution, still vanishes for c = 0): 

kT 

with ©c = , and k the Boltzmann constant. On the other hand, when U 

qt 

is larger than Ut, the transition is first order, and the critical temperature can 
be evaluated numerically. Figure 2. shows the behavior of Tc vs. U in the 
two regions. The value U = A correspond to the vanishing of both the critical 
temperature and cq. 
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Figure 2. can be compared with the rough estimate of the critical temper- 
ature of the long-range order phase whose existence is proved for the Falicov- 
Kimball model in [8] . If one assumes that the phase with could possibly 
be the long range order phase, the qualitative behavior of Tc is in agreement 
with that given by Kennedy and Lieb for large U, whereas it is in contrast with 
the latter for vanishing U. One should notice however that our approximation 
is expected to be more realistic for finite U. 

Finally, in figure 3. we give the behavior of c vs. T for various fillings, still 
for a symmetric state {n = d). The figure shows that the transition is present at 
different fillings, again in agreement with the features of the long-range ordered 
phase described in [8]. 

4. Conclusions 

In this paper we have proposed an improvement of the Fermi-linearization 
technique for electron systems, based on the requirement that the coefficients of 
the linearized operators are Clifford-like variables, with their anticommutator 
equal to an unknown constant c. As an example, we applied such method 
to the Falicov-Kimball model, also giving the self-consistency equation which 
determines the unknown. The latter turned out to behave as a true order 
parameter, which at T = and at half- filling was shown to coincide with the 
discontinuity in the single particle occupation number at the Fermi level in the 
Gutzwiller approximation to the Hubbard model. The behavior of c was thus 
related to the existence of a metal-insulator transition, which again was shown 
to coincide at T = with that hypotized by Brinkman and Rice'-'^^^ . 

The above results suggest that our approximation could be a natural ex- 
tension of the Gutzwiller approach to the case T ^ 0. They also provide a 
physical interpretation to the method, which consists in replacing the hopping 
term by a term which locally still allows to create and annihilate electrons, 
but with an amplitude proportional to the discontinuity in the single particle 
average number at the Fermi surface. 

Moreover, as opposite to the case in which the coefficients of the linearized 
operators were Grassmann variables, the present approximation produces non- 
trivial results even in the case U = 0. 
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This paper was intended as a presentation of the method, and httle efforts 
were devoted to the numerical results in the various cases. Nevertheless, in view 
of the promising results obtained, work is in progress in order both to provide 
a complete phase space at T = and to discuss the T ^ behavior of the 
physical quantities. We expect that also in this case the use of a cluster Bethe 
version of our approximation^^] should give more accured quantitative results. 

Finally, let us stress that the method is of further generality. In particular, 
we expect that it can be straightforwardly applied to the conventional Hubbard 
model, as well as to generalized Hubbard models which have been proposed for 
the study of high-Tc superconductivity. 

The authors gratefully acknowledge stimulating discussions and a careful reading 
of the manuscript by Mario Rasetti. One of the authors (A.M.) also thanks Dan Mattis 
for letting her know about ref. [4]. 
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Figure captions 

kT 

Fig. 1. c vs at different U values: U — (continuous line), U — 1 (dashed 

line), U = 2 (dotdashed line). 
kT 

Fig. 2. — - vs U: continuous line represents second order transition, dashed line 

first order transition. 

kT 

Fig. 3. c vs — at U = 1 and different fillings, in the neutral case {n = d): n+d = 1 
qt 

(continuous line), n + d — .8 (dashed line), n + d — .6 (dotdashed line). 
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